Jacobi vector fields for Lagrangian systems on algebroids 



Michal Jozwikowski* 
Institute of Mathematics, Polish Academy of Sciences 
Sniadeckich 8, P.O. Box 21, 00-956 Warszawa, Poland 
mj ozwikowskiOgmail . com 

March 4, 2013 



Abstract 

We study the geometric nature of the Jacobi equation. In particular we prove that Jacobi 
vector fields (JVFs) along a solution of the Euler-Lagrange (EL) equations are themselves 
solutions of the EL equations but considered on a non-standard algebroid (different from 
the tangent bundle Lie algebroid) . As a consequence we obtain a simple non-computational 
proof of the relation between the null subspace of the second variation of the action and the 
presence of JVFs (and conjugate points) along an extremal. We work in the framework of 
skew-symmetric algebroids. 
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1 Introduction 

Jacobi vector fields (JVFs) and conjugate points appear naturally in Riemannian geometry in 
connection with geodesies. There are at least three different ways to speak about JVFs. Consider 
a Riemannian manifold (M,g) and let £ : [to^i] — >■ M be a geodesic. 

(A) First of all, we may define a JVF in geometric terms: by means of the Levi-Civita connection 
V and the Riemannian curvature tensor R(-, •) of (M,g). It is a vector field J(t) along £(t) 
which satisfies the following second order differential equation 

V^J + #(J,C)C = 0. (1.1) 

(B) In fact (|l,ip (known as the Jacobi equation) appears naturally when considering the tangent 
prolongation of the geodesic equation V^C = 0. Strictly speaking if t t— > ((t,s), where s G 
(—e,e) and C(^0) = C(*)j 1S a one-parameter family of geodesies, then J(t) := 9 s | s _ £(i, s) 
satisfies (|3.4p . Conversely, every solution of (jl.ip may be constructed in this way. This gives 
the second interpretation of JVFs. 
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(C) Finally, JVFs appear in connection with second variation of the energy functional S(C) = 
\ Ito ^(^> Precisely, if the first variation SS(() vanishes (i.e. £ is a geodesies), then the 

second variation reads as (see |12| ) 

5 2 S(()(X, Y) = - jT 1 g (y, V C V^X + R(X, ()() dt - g (y(t), V^(t)) [\ (1.2) 

where X and V are vector fields along £ (and at least one of them is vanishing at the 
end-points). Inside the integral we recognize the left-hand-side of the Jacobi equation. 

This third interpretation of JVFs becomes more important after introducing the notion of 
conjugate points. We say that C(io) an d C(*i) are conjugated along £ if there exists a non- 
vanishing JVF J(t) along £ which vanishes at to and t± (it corresponds to the existence of a 
family of geodesies emerging from C(^o) an d meeting again in C(^i) ( U P to first order terms)). 
The celebrated results of Morse |12| links the presence of conjugate points with the properties of 
the action S(-). 

Theorem 1.1 (Morse). Consider S 2 S(C)(- , •) as a symmetric bilinear form on the space of vector 
fields along £(i) which vanish at to andt\. Then, 

1. The index (i.e. the dimension of the negative space) of 5 2 S(C){',') * s always finite and 
equals the number of points (counted with their multiplicities) conjugated to C(^o)> which lie 
on ( between C(^o) an d C(^i) - 

2. A vector field X(t) along (^{t) is a null vector of 5 2 S(C){-,-) if and only if X{t) is a JVF 
vanishing at to andt\ (i.e. C(^o) and (,{ti) are conjugated along 0. The dimension of the 
null space equals the multiplicity o/C(io) an d C(^i) as conjugate points. 

In fact, much of the theory of JVFs and conjugate points can be rewritten for an arbitrary 
Lagrangian system L : TM — > M. There is, of course, sense in speaking about a tangent 
prolongation of the Euler-Lagrange (EL) equation, i.e. defining a JVF as the first order variation 
associated with a one-parameter family of solutions as in (B). We may also link such JVFs with 
the properties of the second variations of the action functional S(() = J^ 1 L(()dt as in (C). 
Also the notion of conjugate points can be straightforwardly generalized to this situation but, in 
general, we cannot expect any results of the same type as in the first part of Theorem 11.11 unless 
some additional assumptions are put on L. However, the second part of this theorem, i.e. the 
link between existence of conjugate points and the null subspace of the second variation 5 2 S, can 
be maintained. 

What seems to be unclear is the geometric nature of these generalized JVFs (the description 
corresponding to (A)). In this paper we aim to give an answer to this problem. In fact we 
will consider Lagrangian systems of a more general class, namely Lagrangian systems on skew- 
symmetric algebroids. 

Algebroids have gained much attention in recent years providing a useful geometric language 
to speak about mechanics (see [2] and the references therein). However, some of the studies of 
these objects seem to be just vain straightforward generalizations of the classical concepts. We 
believe that this is not the case of this paper. Our main result states that the Jacobi equation 
associated with a Lagrangian system L : E — > R on a skew-symmetric algebroid r : E —> M 
is strictly connected with the EL equations for the lifted Lagrangian system d-r-^ : TE 1 — > M. 
on the tangent lift algebroid Tr : TE — » TM. Note that even starting form the classical case 
of a Lagrangian system on the tangent bundle Lie algebroid tm ■ TM — > M we end up with 
a Lagrangian system on the algebroid Ttm ■ TTM — > TM which is not the standard tangent 
bundle Lie algebroid (in particular, we can interpret (jl.ip as a projection of the EL equation from 
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Ttm '■ TTM — > TM to tm '■ TM —¥ M. This example shows that skew-symmetric algebroids are 
the proper class of objects which allows to speak about JVFs even if we wish to restrict ourselves 
to system on tangent bundle Lie algebroids only. 

To formulate our main results we need to introduce the notion of an admissible curve 7 : 
[to,ii] — > E on an algebroid r : E — > M and the notion of an infinitesimal admissible variation 
5^7 £ T^E of 7 generated by £ G E TCyy . Precise definitions can be found in Section [2j At 
this point an intuitive interpretation will suffice. Treating E as a, bundle of generalized velocities 
(configurations space) of a physical system we may think of 7 as a true physical trajectory of this 
system. The variation ^7, in turn, can be understood as a variation at the level of generalized 
velocities associated with a virtual displacement in the direction of the generator £. 

Theorem 1.2. Consider a Lagrangian system L : E — >■ R on a skew-symmetric algebroid r : 
E — > M and an admissible curve 7 being a solution of the EL equation for this system. A 
curve 7] € E T0 ^ is a JVF along 7 if and only if the associated infinitesimal admissible variation 
<5 r; 7 G T^E is a solution of the EL equation for the lifted Lagrangian system d^L : TE — > R on 
the tangent lift algebroid Tr : TE TM. 

What is more, from Theorem 11.21 and the standard relation of EL equations with variational 
principles, we almost tautologically obtain the relation between the second variation of the action 
functional and the existence of JVFs and conjugate points. This time we have to restrict ourselves 
to a certain subclass of skew-symmetric algebroids, namely almost-Lie (AL) algebroids. 

Theorem 1.3. Let L : E — s> R be a Lagrangian system on an AL algebroid r : E — > M. Consider 
the action ^ 

7 : [to, h]^E .— ► S L (l) ■= [ 1 L( 7 (t))dt 

Jt 

defined on the space of admissible curves. If ' j is a solution of the EL equation then the second 
variation 5 2 S 1,(7) (-, •) is a well-defined bilinear form on the space of curves € E T0 ^^ x 

E® o ^^ t y Moreover, rj(t) € E TQ1 is a JVF along 7 if and only if 

S 2 S L { 1 ){r 1 ,t)=Q for any £ G E° OJ . 

Moreover, if E is a Lie algebroid, then 

5 2 S L ( 7 )(-,-):^o 7 x^ 07 ^R 

is a symmetric bilinear form. The dimension of the null-space of this form equals the multiplicity 
of points j (to) and j(ti) as conjugated points along 7. 

To sum up, our main purposes in this publication are: 

• Defining the Jacobi equation (and its solutions, JVFs) as the tangent prolongation of the 
EL equation on a skew-symmetric algebroid. 

• Giving a geometric interpretation of the Jacobi equation in terms of the EL equation on 
the tangent lift of the initial algebroid. 

• Studying the relation between second order variations of the action functional and the 
existence of JVFs and conjugate points. 

• Doing all the above in purely geometric terms avoiding working in local coordinates when- 
ever possible. 

We are aware that some readers may prefer coordinate rather than geometric description. For 
them we summarize the most important formulas in Section \5\ 
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2 Geometric preliminaries 



In this section we introduce some basic concepts associated with the notion of a skew-symmetric 
algebroid and with the notion of the tangent lift of a tensor field on a manifold. At the end we 
show how these two constructions join in the concept of the tangent lift of an algebroid. 

Notation. We will use standard notions and conventions of differential geometry. Throughout 
the paper we will be working with a vector bundle r : E — >■ M, its dual tt : E* — > M and related 
objects, such as TE, TE*, TTE, etc. We will use YE = ker{TY : TE ->• TM} to denote the 
subbundle of vectors in TE vertical w.r.t. r : E — >• M. By £ 6 E x or G ^Wt) we wm denote 
a section of along a base path x(t) € M (such curves will always be defined on the interval 
[to, ti\). Without some abuse of notation we will denote by E x the set of all such sections (along 
a fixed base path x(t)) and by E® all such sections vanishing at the end-points (i.e. £(to) an d 
£(£i) are null vectors w.r.t. the vector bundle structure r : E — > M). If a manifold carries 
more than one vector bundle structure (like in the case of TE) we will emphasize which vector 
bundle structure we are speaking about. For a tangent bundle TM we will write rather T X M 
than (TM) . We do not introduce any local coordinates at this stage. Readers interested in the 
coordinate description of the theory should read the text parallel with Section [5] where such a 
description is given. 

2.1 Algebroids 

Definition 2.1. A skew- symmetric algebroid (or shortly a skew algebroid) structure on a vector 
bundle r : E — > M is a skew-symmetric bilinear bracket [-, •] on the space Sec(r) of sections of 
r, together with a vector bundle morphism p : E — >• TM (the anchor map) which satisfy the 
Leibniz rule: 

[XJ-Y] = f[X,Y] + p(X)(f)Y (2.1) 

for every X, Y € Sec(r) and / € C°°(M). If, in addition, p maps the algebroid bracket into the 
standard Lie bracket on TM, i.e. 

p([X,Y]) = [p(X),p(Y)] TM , (2.2) 

we call E an almost-Lie (AL) algebroid. If E is AL and, moreover, 

(Sec(r), [•,•]) is a Lie algebra, (2-3) 

we speak about Lie algebroids. 

Algebroids as DVB morphisms. There are several equivalent definitions of the algebroid 
structure on E. For example, it can be understood as a linear bivector field II on the dual bundle 
7r : E* — » M. Such a II defines a double vector bundle (DVB) morphism LT : T*^* — > TE* given 
by contraction. More information on LT, how it is constructed, its coordinate description and 
how the additional conditions (e.g. E being Lie or AL, etc.) are reflected in it can be found in 
[21 [3]. Let us only remark that II is a linear Poisson tensor if and only if E is a Lie algebroid. 

Composition of LT with the canonical DVB-isomorphism TZe ■ T*E — > T*E* (for more infor- 
mation on TZe consult the Appendix) gives a DVB-morphism e := Ilo7^£; : T*E — > TE* over 
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the identity on E*: 

T*E ^TE* . (2.4) 




This morphism gives another characterization of the algebroid structure on E. It will be 
particularly useful in the context of Lagrangian mechanics. 

Relation k. Since the dual bundles of tte ■ T*E — > E and Tir : TE* — »■ TM are, respectively, 
te ■ TE — > E and Tr : TE — > TM, the dual to the DVB morphism e is a relation n : TE — > TE. 
It is a uniquely defined smooth submanifold k C Ti? x Ti? consisting of pairs (X, Y) such that 
p(T E (Y)) = Tt{X) and 

(X,e(A)) Tr = (Y,A) TE 

for any ^4 G T^y^E 1 , where -) Tt . and are natural pairings. This relation can be put 

into the following diagram of "double vector bundle relations": 

TE t TE . (2.5) 




Canonical DVB morphisms kq and olq Of particular importance is the relation k associ- 
ated with the canonical tangent bundle Lie algebroid structure on TQ. We will denote it by kq 
to emphasize where form it originates. In fact, kq is not just a relation but a DVB-isomorphism 
which intertwines the two canonical projections ttq : TTQ — >■ TQ and Ttq : TTQ — > TQ: 



TTQ ^ TTQ 




A nice geometric interpretation of kq is the following. A vector A € TyTQ can be represented 
by a curve s\ i-> Y(s\) £ T^^Q with 1^(0) = Y. Since every vector itself Y(si) can be 
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represented by a curve S2 q(s\, S2) G Q where q(s\, 0) = q(s\), the vector A can be regarded 
as an infinitesimal object associated with a homotopy (s\,S2) i-> q(s\,S2) G Q- The projection 
Y = t tq(-<4) € TQ is an infinitesimal object associated with the curve S2 1— >■ <7(0, S2), while the 
projection X := Ttq(A) G TQ with the curve si h-> g(si,0). Now we can reverse the order of 
variables si and S2 to obtain a new homotopy g(si,S2) = q(s2,si). The vector # corresponding 
to q is the image of A under kq. Obviously, ttq(A) = Ttq(B) and Ttq(A) = ttq(B). Moreover, 
we see that kq is an involution (hence isomorphism). 

Relation kq is dual to a DVB morphism eq : T*TQ — > TT*Q. Since kq is an isomorphism, 
so is eq. Its inverse is usually denoted by ckq : TT*Q — > T*TQ: 



T*TQ ■< TT*Q 




It is defined via the equality 

(X,y) TTQ = (K Q (X),a Q (V)) TTQ , (2.6) 
for X G TTQ and ^ G TT*Q such that the first parring makes sense. 



Admissible curves and infinitesimal admissible variations. If r : E — > M carries a skew- 
algebroid structure, then thanks to k we may introduce the notion of infinitesimal admissible 
variations of curves in E. We begin with the definition of an admissible curve. 

Definition 2.2. A curve 7 : [to, t\] — > E is called admissible if the following condition is satisfied 

= 4(To 7 (t)), far every t€ Mi], (2.7) 
at 

i.e. the tangent prolongation of the base projection of 7 agrees with the anchor of 7. 

In the standard case E = TM admissibility of 7 means that 7 (t) = C(0> wh ere C := r 7, that 
is 7 is a tangent prolongation of its base projection 

Let now 7 : [to,ii] — > E be an admissible curve. A curve £7 : [to, ^l] ~^ TE which projects to 
7 under te '■ TE — >■ E is called an infinitesimal variation 0/7. Among all variations of a given 7 
we distinguish the subclass of infinitesimal admissible variations. These are constructed in the 
following way. 

Take any path £ : [to, ti] — > £7 which projects to the same base path as 7, i.e. r o £ = r o 7. 
The tangent prolongation £(t) is a vector field in TE along £. Now there exist an unique vector 
field ^7 in Ti? along 7 which is K-related to £(t), that is 

<% 7 (t) = k 7 (e(t)) . 

The lower index k 7 indicates that ^7 G T 7 £\ Note that £ 1— > <5^ 7 is linear and that from (|2.5p 
it follows that ^7 projects to p(£) under Tr: 

Tr(5 a ) = p(£). 

We can summarize the above construction in the following definition: 
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Definition 2.3. Infinitesimal variations of the form ^7 are called admissible. The curve £ is 
called an infinitesimal generator of We distinguish admissible infinitesimal variations with 
vanishing end-points generated by £'s such that G E T0 ^r t A for i = 1,2 is a null vector in 

t:E^M. 

The intuitive interpretation of admissible curves and admissible variations is the following. 
We may treat E as a bundle of generalized velocities over M with true velocities on M obtained 
after applying the anchor map p : E — s> TM. Admissible curves are these curves in the generalized 
configuration space E which correspond to true movements. A generator £ : [to^l] E of an 
admissible variation can be now understood as a generalized direction in which a variation is 
performed and the full infinitesimal variation 5^ is a variation in the generalized configuration 
space corresponding to the movement in the direction of a generator. 

Remark 2.4. It may be useful for the reader to describe in detail the connection between relation 
k and infinitesimal variations in the particular case of k = Kn. Let namely t 1— > q(t,s) be a s- 
parameterized family of curves in Q. It generates a variation Y(t) = d s \ Q q(t,s) G T q ^Q along 
q(t) := q(t,0). Denote its i-derivative by A(t) := dtY(t) G TymTQ. On the other hand, we 
can lift q(t,s) to curves X(t,s) = dtq(t,s) G T q ( t ^Q. Again t *— > X(t,s) is an s-parameterized 
family of velocities. It defines a variation B(t) = d s \ Q X(t, s) G Ty(t)TQ of a velocity curve 
X(t) := X(t,0) = d t q(t) G T q ^Q. Now B(t) ^ A(t), as the first projects to X(t) and the other 
to Y(t) under ttq. We have, however, KQ(A(t)) = B(t). 

In other words, KQ(dtY(t)) is a variation at the level of velocities associated with a variation 
Y(t) at the level of base trajectories. 

As we have seen above, on the tangent bundle Lie algebroid infinitesimal admissible variations 
are infinitesimal objects associated with true variations, i.e. one-parameter families of admissible 
curves. It turns out that this correspondence holds if E is an almost Lie algebroid (in particular, 
it holds for every Lie algebroid) . It was first observed in [2] (see also [1] for a detailed discussion) . 

Lemma 2.5. Let r : E — > M be an AL algebroid. Then the tangent space of a Banach manifold 
ADJ\A([to,ti],E) of admissible paths 7 : [to?^l] ~~ ► E equals 

T 7 ADM([toM E ) = {^7 : £(*) G E TOl(t) }. 

Let us note that this nice interpretation is not true for general (not AL) skew-algebroids. In 
fact almost-Lie algebroids are precisely these skew-symmetric algebroids which posses the above 
property (see the references above). 

2.2 Tangent lifts 

Now we will introduce the notion of a tangent lift of a geometrical object on a manifold and 
apply it to a skew-algebroid. 

Tangent lifts of multivector fields A tangent lift dx of a multivector field on a manifold 
Q was introduced in [6] as a natural extension of a derivation die of differential forms (in fact, 
it can be extended to the whole tensor algebra T°(Q)). Roughly speaking, die is a differential 
operator which maps tensor fields of type a (a is a multi-index) on Q to tensor fields of the same 
type on TQ. 

d T : T [a \Q) ^T (Q) (TQ). 

For our purposes it will be enough to define dx on multivector fields. We will state only 
these results which are important to our considerations, referring to [6] as an extensive source of 
information. 
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Introduce local coordinates (q a ) on Q and induced coordinates (q a ,q^) on TQ. For a multi- 
vector field X G 3T(Q) where X(q) = X ai - ar (q)d q *i A ... A d q <* r , its lift d T ^ is locally defined 
by 

d T X{q, q) = —-}*Lqf>d r i A . . . A + ^ X ai - ar (q)d#*i A ... A d q a m A ... A d q « r (2.8) 

m=l 

The general geometric definition can be found in [5]. For us it will be sufficient to give it just in 
cases of fields of order r = 0, 1 and 2. 

Let / G C°°(Q) be a smooth function (a 0- vector field). We have simply 



d T /(X) = (df(q), X), where X G T g Q; 



(2.9) 



that is, d-r/ is a linear function on TQ defined by the differential d/. 

A vector field V G 3t(Q) can be understood as a map V : Q — > TQ covering idg. The map 
(ItV corresponding to the lifted vector field d^V is a composition TV o kq : TQ — s> TTQ 



TQ 

id 

TQ 



TV 



kl 



^ TTQ 

TrQ , 
— TQ 



id 



^ TTQ 
-+TQ 



(2.10) 



Another way of thinking about d^V is the following. Let <&* : Q — > Q be the (local) flow of 
V. The map T<E>' : TQ — > TQ is also a (local) flow, hence it corresponds to some vector field on 
TQ. This field is precisely dx^- 

Finally, a bi-vector field IT € X 2 (Q) can be associated with a vector bundle morphism II : 
T*Q — > TQ given by contraction. The map diTI : T*TQ — > TTQ corresponding to dxll G 
j£ 2 (TQ) is defined by means of the following diagram: 



T*TQ 
TT*Q 



d T n 



Tn 



TTQ . 

KQ 

'Q 



(2.11) 



2.3 Tangent lift of an algebroid 

As was mentioned before, a skew- algebroid structure on r : E — > M is a linear bivector field 
n G X 2 (E*). Now the tangent lift d T II is a bivector field on TE* linear w.r.t. Tvr : TE* ->■ TM. 
Consequently, it defines a skew-algebroid structure on Tr : TE — > TM. We will call it a 
tangent lift of an algebroid on E and will denote it by d^E. The reader can consult [5] for more 
information on d^E. 

As has been shown in [6], dT commutes with the Schouten bracket, hence if IT was Poisson 
then so is diTI. Consequently, a tangent lift of a Lie algebroid is again a Lie algebroid. One can 
show that the lift of an AL algebroid is an AL algebroid itself. 

In particular, TTM posses two different Lie algebroid structures: the standard tangent bundle 
structure on t tm : TTM -> TM and the lifted structure on Tt m : TTM -> TM. The two are 
isomorphic via the map km- 
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The structure of d-^E. We shall now describe the structure (the morphism e, the anchor 
map, the relation n) of the tangent lift Lie algebroid d^E. We will denote these objects by d^e, 
dx/f and d^K, respectively. 

Theorem 2.6 (0). For the tangent lift algebroid d^E the morphism d^e : T*TE — > TTE* is 
defined by means of the following diagram 



T*TE ■ 



TT*E 



Te 



TTE* 



TTE* 



(2.12) 



The anchor map d^p : TE 
TTE — > TTE is given by 



TTM is the composition dxp = k-m ° Tp and the relation d^K '■ 



T 



E- 



-T r 



TTE 



Tk 



E . 



(2.13) 



TTE 



Proof. To prove (|2.12p observe that, by definition, d-^e is the composition of TZte '■ T*TE — > 
T*TE* with d^A : T*TE* — > TT.E*. Composing two commutative diagrams we get the follow- 
ing: 

Kte mm „ <1tA 



T*TE ■ 



a E 



TT*E- 



E* 



TT*^?* 



TA 



TTE* 



TTE* 



Since the bottom row of the above diagram is TA o TTZe = T(A o TZe) = Te, we have (|2. 12[) . 
The formula for d^p is easily obtained by projecting d-pe to the base: 



Te 



TT*£ — 



T*TT E 

TE=— 



OLE 



id 



■ T*TE 



TT TE 

TE — 



d^e 



d T p 



TTE* 



'<E 



TTtt 

TTM 



TTE* 

TTtt 

: TTM 



Finally, ()2.13p can be proved as follows. By definition, d^n is the relation dual to the 
DVB morphism dx£- On the other hand, one can show that the pair of isomorphisms (ke,k e ) 



preserves the parring 



'TTr 



: TTE x T TM TTE* 



and that, by definition, the pair of 



isomorphisms (k e ,a E ) maps the parring (■,•) : TTE Xte T*TE — > R to the parring 



<•, -) T : TTE x TE TT*E -> R (equation §ZG§ ). It follows that relation Kg 1 o dxft ° k-e is dual 



to the morphism k e , o dp£ ° = Te. On the other hand relation Tk is dual to Te, hence 
Tk = k^ 1 o d^K o ke, which is precisely (I2.13p . □ 



3 EL equations and JVFs 

Now we will recall (after [2j [3]) the geometric construction of the EL equations for a Lagrangian 
system L : E — > R on an algebroid r : E — > M. Later we will define JVFs for such a system as 
the solutions of the tangent lift of these equations. Using the machinery introduced in Section 
[2j we will be able to interpret JVFs as solutions of the lifted Lagrangian system dx^ : TE — > R 
on the tangent lift algebroid Tr : TE — >• TM. 
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Lagrangian dynamics By a Lagrangian system on an algebroid r : E — > M we will 
understand a smooth function L : E — > R. Given such L we can construct two maps 
\l = te* ° e ° dL = T*r o dL : E — > E* (the Legendre mapping) and = e o dL : E -> TE* 
(the Tulczyjew differential) as shown on the diagram 

T*E - ^TE* . (3.1) 



dZ 



E---- L --^E 



We define the Euler- Lagrange {EL) equation for a curve 7 : [to, t\[ — >• E as follows: 

A(A L ( 7 (t))) = A L ( 7 (t)). (3.2) 

In other words 7 (i) is a solution of the EL equation if A^( 7 (t)) is the tangent prolongation of 
Ax(7(*))- 

Observe that every such 7 is automatically admissible. Indeed, we need a simple diagram 
chasing on f|2.4() to see that T7r o A^( 7 ) = p( 7 ) and 7r o A £,(7) = t o 7 , hence on the base level we 
get 

4 (r o 7 (t)) = A (vr o A L ( 7 (t))) = o ( ±\ M t))) <P Tvr o A L ( 7 (t)) = p( 7 (*))- 



dt dt v v v /7/ V dt 



Note also that for any admissible curve 7 (not necessarily a solution of (|3.2p ) the difference 
lt< 



A^( 7 (t)) — ^(Ax,( 7 (t))) is a vertical vector. Indeed 



Tvro (a l ( 7 (*)) - ^A L ( 7 (t))) = Tvro A L ( 7 (t)) - ^ (tt o A L ( 7 (t))) = p( 7 (*)) - ^ (r o 7 (i)) = 0. 

Variational interpretation of EL equations. Consider a natural action functional <Sl( ) 
associated with a Lagrangian function L : E — >• R: 

7 : [to, *i] E 1 ^ Sl(t) := / 1 L( 7 (t))dt. 

•no 

The first variation of Sl(') at 7 in the direction of an infinitesimal variation £7 is defined as 
(dS L (7),<5 7 ) := f h (dL( 1 (t)),5 1 (t))dt= f * d T L(&y(t))dt. 

J to J to 

The following result is according to Grabowska and Grabowski [2|. 

Theorem 3.1. For an admissible curve 7 : [to,ti] — ^ E an d every infinitesimal admissible 
variation ^7, the following equality holds: 



to 



(d&( 7 ),«€7> = Z* 1 <e(i),v Ai(7(t)) _d AL(7W) >dt + (e(t),A L ( 7 (t)); 



(3.3) 

to 



//ere u\\_A_\ L OyM) S E* denotes the vertical part of the vertical vector A/, (7(f)) 

dt' 



d rAz,(7(i)) ^ VE* C TE*, i.e. its image under canonical identification YE* ~ E*. 



Remark 3.2. Observe that fj3.3[) implies that for an admissible curve 7 : [to, ^l] ~~ * E the 
following conditions are equivalent: 
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1. The curve 7 is a solution of the Euler-Lagrange equation for the Lagrangian system L : 
E^Kon the skew-algebroid r : E — >• M, 

2. For every infinitesimal admissible variation ^7, we have 

(ds L ( 1 ),6^) = (at),x L h(t))) 

to 

3. For every infinitesimal admissible variation ^7 with vanishing end-points, we have 

(dS L ( 1 ),5 a )=0. 

Remark 3.3. The above result can be understood as a variational formulation of the EL equa- 
tions. Note, however, that it is in some sense different form the "natural" variational formulation. 
In the "natural" approach an admissible curve ^(t) is an extremal of the action Sl(-) if for every 
1-parameter family (homotopy, variation) of admissible curves 1 1— > ^(t, s) such that j(t, 0) = 7(f) 
(and perhaps satisfying additional conditions such as fixing the end-points) a point s = is a 
critical value of the function s 1— > <Sl(7(-, s)), i.e. 

(dS L ( 1 ),d s \ Ql (; S ))=0. 

In Theorem 13.11 and Remark 13.21 we present however the "formal" approach - we assume 
vanishing of dSi,^) on infinitesimal variations of 7 of certain type (^7 - infinitesimal admis- 
sible variations with vanishing end-points) without assuring ourselves that these variations are 
generated by homotopies. Of course from Lemma [2. 5 1 we know that if E is an AL algebroid then 
every infinitesimal admissible variation is generated by a homotopy, and hence the "natural" and 
the "formal" approach coincide. In general (if E is not AL) the geometric EL equations (|3.2p are 
only "formally" variational. 

JVFs. Now we will give the most important definition of this paper. To motivate it, consider 
an s-parameterized family of solutions of (|3.2p 1 1— > j(t, s). Then, 



_d 

d^ 



^A L ( 7 (i, s)) = ^A L ( 7 (i, s)) = TA L (d s \ 7(t, «)). 



Now, after using ke* to inverse the order of derivatives on the left-hand side, we get 



«^oTA L (a.| o7 (*, S )) = ~ 



Al(7(M)) = £(TA l (0 s | o7 (M))) 



dt 



that is J(t) := d s \ r y(t,s) satisfies 

^TX L (J(t)) = K E * o TAx(J(t)). (3.4) 

In other words, the above equation is the tangent prolongation of (|3.2p . Obviously, it is an 
implicit linear differential equation of the first order. It is tempting to call (|3.4p Jacobi equation, 
yet this notion would not agree with (jl.ip . The reason is that in the classical geodesic problem 
we treat the geodesic equation not as the equation for an admissible curve 7 = £ 6 TM but as 
an equation for the base curve £ G M. Therefore we do not look for arbitrary variations of 7 = £ 
but only for admissible ones (i.e. variations generated by the true variations on the base). This 
leads us to the following formulation. 
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Definition 3.4. Let 7 (f) be a solution of the EL equation fj3.2j) for a given Lagrangian system 
L : — > K on a skew-symmetric algebroid £\ A Jacobi vector field (JVF) along j(t) is a curve 
7?(t) € E rol n\ such that the associated admissible variation 5 v j satisfies (|3,4p . i.e. 

^TA L (^ w7 (i)) = ke* o TA L (<f„ (t)7 (i)). (3.5) 

Equation (|3.5p will be called the Jacobi equation. 

We say that the points 7^0) and 7 (ti) are conjugated along 7 if there exists a non-zero JVF 
rj(t) € E TO jM such that 77(^0) = vi^i) = 0. Obviously, such JVFs form a vector space. Its 
dimension will be called the multiplicity of conjugate points 7 (to) an d 7(ti). 

Observe that, since r/ 1— > <5^ 7 is linear and (|3.4p is linear in J, (|3.5p is a second order linear 
implicit differential equation for 77. To see that (|3,5p coincides with in the Riemannian case 
see Section [6J 

Remark 3.5. Let us now discuss the admissibility of solutions of (|3.4p and (|3.5|) . Since every 
solution of (|3.2p was admissible, every solution J 6 T 7 i? of its tangent prolongation (|3.4p belongs 
to T^AT>A4([tQ, t\], E). Indeed, an application of TT-7T to (|3.4p gives 

j t (Tr o J(t)) = km o Tp(J(T)) = d T p(J(t)), (3.6) 

which is the tangent prolongation of (|2,7|) and hence the condition for J € T 7 .E to belong to 
T 7 AOM([i ,ii],#). 

Now, if £7 is AL, then (Lemma 12. 5p every solution of (|3.4|) is of the form J = 5 v j. In general, 
if E is not AL, then there can be no solutions of (|3.4p of the form J = 5^ 7 . Indeed, recall that 
Tr o 5^7 = p(n), hence ^ (Tr o 5^7) = Tp(rj). Consequently, 

Tp(r)) = (Tr o ^7) = k m ° Tp(5r,j) = k m o Tp(K 7 (r))) 

(i.e. J = 5^7 satisfies (|3.6p ) if and only if 

k m oTp(r)) = Tpo K(i]). 

The condition km Tp = Tp o k means that Tp maps k to km, that is p is a morphism of 
algebroids E and TM (in other words, -E is AL). 

The main result. It turns out that an admissible variation associated with a JVF is again a 
solution of the EL equation. 



Theorem 3.6. Equation (|3.4p is an Euler- Lagrange equation corresponding to the Lagrangian 
d^L : TE — > M and the tangent lift algebroid structure d^E. Consequently, if rj is a JVF along 
7; then the associated admissible variation ^7 is a solution of this EL equation. 

Proof. Euler Lagrange equations can be constructed via a well-defined geometric procedure as 
described on diagram (|3.ip . Since, by Theorem 12.61 we know how to describe the structure of 
the tangent lift algebroid d^E, the geometric construction of the Euler-Lagrange dynamics for 
d^L on d<rE is the following: 



dd T L 
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On the other hand, fj3.4|) can be described as follows: 



T*TE- 

OiE 

TT*E 



Te 



TTE* 



TTE* 



TdL 



TE 



TXr 



■TE* 



To prove the assertion it is enough to check that a^o (TdL) = d(dxL) and TA^ = Xd T L- To 
justify the first equality take A € TTE represented by a homotopy j(t,s) G E. It is easy to see 
that 



(d(d T L),A) 



_d 

TTE ~ ds 



d 

,di 



L( 7 (t,s)). 



Now, 



[a E oTdL,A) 



TTE 



"9 <Td W .4» TrE = i 



d_ 

dl 



since ke{A) is represented by j(s,t). As the derivatives commute, 



(d(d T L),A) 



= (a E oTdL,A) . 

tte x 1 tte 



Finally, observe that 



TA L = T(t £ .o £ o dL) = Tt e * ° Te o TdL = (Tt e * o k e *) o (k b . o Te) o TdL 

tte* ° dx£ o oiE o TdL = ttb* ° dx£ ° ddxL = Ad T L- 



□ 



4 Second order optimality conditions 

In this section we shall study the relation of JVFs and variational principles. Our main result 
is Theorem 14.71 which shows the link between the presence of conjugate points and null vectors 
of the bilinear form 5 2 Sl - the second variation of the action Sl- The result is proved almost 
tautologically thanks to the interpretation of JVFs and Jacobi equation in terms of EL-equation 
given in Theorem 13.61 Much attention is, however, needed to give a proper definition of the 
second variation of Sl in the case of a Lagrangian system on an algebroid. Therefore we will 
start form general considerations about the second variation of a function at its critical point (or 
a critical point relative to a VB morphism). 

Throughout this section we will be working in the AL algebroid setting. It turns out that this 
is needed for the definition of the second variation to work well. Otherwise we have no control 
of the vectors tangent to the space of admissible curves in E. Moreover, the symmetry of 5 2 Sl 
requires the Lie algebroid setting (see Lemma l4.6p . 



What is the second variation of a function? Consider a smooth function F : Q — > R 
defined on Q (it can be either a finite-dimensional or a Banach manifold). Naively one defines 
the second variation of F at q in the directions T g Q as 



5 2 F(q)(X,Y) 



d 

ds\ 



d 

ds 2 



F{q(si,s 2 )), 
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where q(si,S2) is such that X G T q Q is a vector represented by the curve s% h-> (7(si,0) and 
1" G T g Q is represented by S2 h-> q(0,S2)- In fact, ^ | ^j | -F(g(si, $2)) = d-rd^F^), where 
^4 G TTQ is a vector corresponding to the homotopy (si, S2) 1— > <z(si, ^2)- Observe that since ^ 

and gj- commute, then dTdT-P(*4) = dTdi-P^Q^))- 

To make sure that the second variation of F at q is well-defined, we should therefore show 
that given any other X G TTQ such that t T q(X) = t T q(A) = Y and Ttq{X) = Ttq(A) = X 
we have 

d T d T F(X) = d T d T F(A), 
i.e. that the value of dicd'YF(A) depends on TQ-projections of A only. But, 

d T d T F(^4) - d T dT-F(^') = (dd T F, A-X), 

where we subtract w.r.t. the vector bundle structure ttq on TTQ. Since A and X have the same 
T TQ-projections, A — A G TyTQ is vertical w.r.t. Ttq, and hence in the canonical identification 
VTQ ~ TQ Xq TQ it corresponds to a pair (Y, Z). A simple calculation shows that 

(dd T F,A-X) = d T F{Z) = (dF{q),Z). 

Consequently, 5 2 F(q) is well-defined if and only if dF(q) = 0, i.e. q is a critical point of F. 
We can summarize our considerations as follows: 

Proposition 4.1. Assume that q G Q is a critical point of F. Then the second variation of F 
at q in the directions X, Y G T g Q is well-defined via the formula: 

5 2 F(q)(X,Y):=d T d T F(A), (4.1) 

where A G TT g Q is any vector projecting to X and Y under Ttq and ttq, respectively. 
Moreover the second variation is symmetric, i.e. 

5 2 F(q)(X,Y)=5 2 F(q)(Y,X). 

The second variation at a relative critical point. From the point of view of Variational 
Calculus it is more natural to consider critical points of a function relative to a vector bundle 
morphism, rather than just critical points. Our geometric setting is now the following: F : Q -4 R 
is a smooth function, a : S — > N is a vector bundle and, additionally, we have a vector bundle 
morphism P : TQ — > S over p : Q — > N. A point q G Q is called a critical point of F relative to 
P (which will be denoted by q G Critp(-F)) if 

(dF(q) , X) = for every X G T q Q n ker P. 

To give a concrete example think about the standard variational problem associated with a 
Lagrangian L : TM — )• R. We may take Q = ,4Z?.A4([to) ti], TM) to be the space of admissible 
curves in TM (i.e. tangent lifts of smooth curves ( : [t ,*i] M); N = M x M , a = t n : TN -4 
N, and P = Tp : TQ -> TiV, where p sends a curve C(-) to its end-points (C(*o), C(*l)) G M X M. 
Define F(C) := L(C(t))dt. Now the solut ions of the EL-equation (extremals of the action) are 
precisely critical points of F relative to Tp (dF vanishes on infinitesimal admissible variations 
with vanishing end-points - cf. Remark 13.21 and Theorem 13.11) . 

Now we would like to define the second variation of F at q G Critp(i ? ). Obviously, the 
idea from Proposition 14.11 will not work directly. Indeed, taking vectors A,A& TTQ as in the 
previous paragraph we have 

d T d T F{A) - d T d T F(X) = (dF(q),Z), 
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which is non-zero unless Z € kerP. To guarantee this we cannot take arbitrary vectors A and 
A with the same TQ-projections X and Y, but have to restrict ourselves to a certain subclass 
of such vectors. 

It is easy to check that if A — A € TyTQ is a vertical vector corresponding to (Y, Z), then 
TP(_4) and TP(A ) project to P(Y) and Tp(X) under ts and Tcr, respectively. The difference 
TP(.A) - TP(A') G T P(y) TiV is a vertical vector w.r.t. Tcr corresponding to (P(Y),P{Z)). 
It follows that Z G kerP if and only if TP (A) - TP(A') corresponds to (P(Y),0). This in 
turn implies TP (A) = TP(A ). In other words, formula (|4.ip will give a good definition of the 
second variation of F at q G Critp(P) if we consider A's with fixed TP-projections. Of course 
TP(^t) € TS should be consistent with TQ-projections of A, i.e. it should project to P(Y) 
under t$ and to Tp(X) under Tcr. Let us summarize: 

Proposition 4.2. Assume that q G Critp(P). Consider an arbitrary map assigning to a pair 
of vectors X G T p ^N and Y G Sp(g) a vector h(X,Y) G TS such that t$ ° h(X ,Y) = Y and 
To~ oh(X ,Y) = X (i.e. h(X ,Y) is a lift of X to T y S). Then we may define the second variation 
of F at q w.r.t. h via the formula: 

6 2 F(q,h)(X,Y) :=d T d T F(A), 

where A G TTQ is an arbitrary vector projecting to X under Ttq, to Y under ttq and to 
h(Tp(X),P(Y)) under TP. 

In general there is no canonical choice of h unless a is provided with some additional structure. 
For example given a connection on a, we may choose h(X, Y) to be the horizontal lift of the 
base vector X to T y S. Observe that if the connection is linear, then 5 2 F(q, h)(X, Y) becomes a 
bilinear form. Indeed, to prove linearity w.r.t. X observe that if A € TTQ projects to X, Y and 
h(Tp(X),P(Y)); and A' to X' , Y and h(Tp(X'),P(Y)), then the linear combination XA + fiA' 
(taken w.r.t. t T q) projects to AX + /j,X' , Y and Xh(Tp(X), P(Y)) + fih(Tp(X'),P(Y)) = 
h(Tp(XX + fiX'),P(Y)). Similarly, using the VB structure Ttq on TTQ and linearity of h(X, Y) 
w.r.t. Y (the connection is linear) one can prove linearity of 5 2 F(q, h)(X, Y) w.r.t. Y. 

The symmetry of 5 2 F(q, h) may be quite problematic. In general one requires that if A € 
TTQ is a vector projecting to X, Y and h(Tp(X), P(Y)) then kq(A) projects to Y, X and 
h(Tp(Y), P(X)) (i.e. kq preserved the distinguished class of vectors in TTQ). If this is the case 
then the symmetry of the second variation follows from the KQ-invariance of dTdxP: 

5 2 F(q,h)(X,Y) = d T d T P(^) = d T d T F(K Q (A)) = 5 2 F{q,h)(Y,X). 

Natural examples of such a situation appear when considering symmetric connections. Take a 
tangent bundle a = tn : TN — > N; p : Q — s> N to be any smooth map and P = Tp : TQ — > TN. 
If h is a horizontal lift of a symmetric linear connection on tn then kn ° h(X , Y) = h(Y , X) for 
any X, Y G T X N. Now if A G TTQ projects to X, Y and h(Tp(X),Tp(Y)) then kq(A) projects 
to Y , X and 

TP(k q (A)) = TT p (k q (A)) = K N (TTp(A)) = K N (h(Tp(X),Tp(Y))) = h(Tp(Y),Tp(X)), 

and hence 5 2 F(q, h) is a symmetric bilinear form. More generally, one can check that if a : S — >• N 
carries an algebroid structure; P : TQ — > S is an algebroid morphism and h is a horizontal lift of 
a symmetric linear connection on a (symmetric w.r.t. to the algebroid structure), then 5 2 F(q, h) 
is a symmetric bilinear form. 

Without the presence of a connection on a we may define the second variation 5 2 F(q)(X, Y) 
after restricting our attention to Y's belonging to ker P. Since P(Y) = we do not have to choose 
h(X ,Y) for any X and Y, but it suffices to define h(X,0) only. There is the canonical choice 
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h(X, 0) = 8jt - the null vector w.r.t. Ta at X G TA. Since r^-addition gives 9jt + &—' = 8— —> 
and Ta-addition gives + = 6^, the second variation will be a bilinear form for this choice 
of h. 

Of course in this setting there is no sense speaking about the symmetry of the second variation 
unless we restrict ourselves to X, Y G T q Q n ker P. It seems that in general the second variation 
will not be symmetric, without some additional assumptions on P and a. To sum up: 

Lemma 4.3. Let q G Q be a critical point of F relative to P : TQ — > S. Then for X G T q Q and 
Y G T q Q H ker P we may define the second variation of F at q in the directions X, Y by 

5 2 F(q)(X,Y) := d T d T F(A), 

where A G TTQ is an arbitrary vector projecting to X under Ttq, to Y under ttq and to #Tp(X) 
under TP. 

Moreover, 5 2 i ? (g)(-, •) is a bilinear form. 

The second variation in variational problems. We would like to apply the construction 
from Lemma 14.31 to variational problems on AL algebroids. In this case we can consider Q = 
AT>A4([tQ, t\], E) - the Banach manifold of admissible paths 7 : [io^i] — > E and F = Sl to be 
the action associated to a Lagrangian system L : E — > M 

S L (l) ■= L( 7 (t))dt. 
J t 

Since E is AL, the tangent space of Q at 7 consists of admissible variations T 7 Q = {<%7 : 
£ G E TO y}. Assume that £ 1— > 5^ is an isomorphism of vector spaces E rol and T 7 Q. (In general 
£ 1 — y (5^7 is just a linear epimorphism. In such little more attention is needed to construct 

the second variation of Sl - see Remark 14. 5p . An admissible curve 7 is a solution of the EL 
equation <j3.2|) if and only if 7 is a critical point of Sl relative to the map P : TQ — >■ E x E, 
which sends to (£(to)>£(ii))- Indeed, T 7 Q Piker P is precisely the set of admissible variations 
with vanishing end-points. 

Our goal is to define 5 2 Sl (7) (^7, ^7) - the second variation at 7 (which is a solution of 
the EL equation) in the directions of 5 v ^f and 5^. For simplicity we will use the abbreviated 
notation 5 2 Sl{^){ti, £) and will treat the second variation as a function of generators r\ and £. 

To proceed within the scheme sketched in the previous paragraph (in Lemma I4.3[) we need 
more information about the iterated tangent space TTQ. To describe T^TQ consider an s- 
parameterized family of admissible variations ^(t, s )7(t, s) around S^u\j(t) = ^(t,o)7(^ 0). Denote 
by A the differentiation at s = 0. Now 

A(5 €(M)7 (t, s)) = A (k 7M (£(M))) = Tk Aj (A(0) = T KA7 ( KB (A£)-) = 
= K E O («£ o Tka 7 o K E ) ((A£)') = K E d T KA 7 ((A£)') , 

where A 7 (t) = £| 7(M) G T 7 £ and A£(i) = J;| £(i,s) G T ?J B. Observe that A7 and A£ 
project to the same path Ax(t) := Tt o A£(i) = Tr o A^ft) = ~§^\ x{t,s) and that A7 is 
dx-E-admissible, i.e. 

(Tr o A7)' = (Ax)' = k m ° A(x) = k m o A(p( 7 )) = k a/ o Tp(A~f) = d T ^(A7). 

If follows that dx^A7 (A£)' is a dT-E-admissible variation along a di^-admissible path A7 with 
a generator A£. With some abuse of notation we will denote it by Sa^Aj. Since we deal with an 
AL algebroid and j(t, s) was a family of admissible paths, from Lemma 12.51 we have A7 = 5 v ^f 
and Ax = p(r]) for some r\ G E rory . 
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We see that elements of T$^TQ are generated by elements A£ G T^P and rj G E to ^ satisfying 
the compatibility condition Tr o A£ = p{rf). Such an element will be denoted by <5(A£,r/)7- The 
precise formula is 

<5 (A$, J? )7 : = k>e o <5a^(^7)- 

By construction <5(A£,r?)7 projects to 0^7 under tte an d to (^7 under Tte- Let us note that since 
P(<% 7 ) = (£(to),£(ti)), then TP ((^,,7)7) = (A£(i ), A£(*i)). Due to Lemma IP we may define 
the second variation of Sl at 7 after restricting our attention to tte <5(A£,r?)7 = ^§7 £ kerP, 
i.e. £ G P° D7 and TP (^^7) G kerTP i.e. A£(i ) and Af(ti) being null vectors w.r.t. 
Tr : TP — > TM. To sum up we can state the following lemma: 

Lemma 4.4. Let 7 be a solution of the EL equation (|3.2p . Formula 

S 2 S L h)(v,0 ■= d T 2 S L (5 {A ^) = f * (dd T L(^7),*(A€^7(t)>dt; (4.2) 

■/ to 

where A£ G T^P is any curve such that Tr o A£ = Tr o A7 = p(rj) and A£(to), A£(ii) are rm/Z 

vectors w.r.t. Tr : TP — > TM; defines a bilinear form 5 2 Si J ('y)(-, •) : P1-07 x P^o7 — ^- ^ e ^ZZ 
ca// it the second variation of Sl at 7. 
Moreover the following equality holds 

5 2 S L (~,)(ri,0 := f <Ae(t),V K ^ oTAL(5j;7W) _, TAi(5r;7(t)) > Tr dt. (4.3) 

J to 

Proof. The fact that the second variation is well-defined by formula (|4.2p follows from Lemma 
14.31 We need to check that equation (|4.3p holds. First of all using the K^-invariance of dT 2 <Sx 
we get 

8 2 S L (j)(rj,£) = d T d T S L (S^^) = d T d T Sx («£ ^(^7)) = d T d T .Sx (^(^7)) • 

Since dxSx = Sd T L, the later equals dT<Sd T L (5a^(^»?7))- We can use (|3.3p for a Lagrangian 
system dTp : TP ->Ron dxP to get 



<5 2 Sl(7)(",£) =dT^d T i(^A7) = (dS^Lidr,-/), 5 A ^(5 v j)) 



to 



( A ^)> V^M^W^W))^ + ( A ^)> Ad T L(^7(i))> Tr 



ti 



fo 



The last summand vanishes as A£(t) vanishes at the end-points. Finally note that from the proof 
of Theorem 13.61 we know that Ad T z, = ke* TA^ and that Ad T L = TA^. This finishes the proof 
of (|P). □ 



Formula (|4.3p generalizes the classical formula (|1.2p - it relates the second variation of Sl to 
Jacobi equation (|3.5p . The details are provided in Section 16.11 



Remark 4.5. Note that the variational interpretation of the EL equation on algebroids (Remark 
I3.2p is slightly more general than what we considered in the paragraph about relative critical 
points. Solutions of the EL equation are critical points of the map 

MAV{[hM],E) — > TQ ^ R 

£ 1 — > 5^7 1 — > dP(5^7) 

relative to the map P : MAV([t , ti], E) 3 £ 1 — > (£(*o), G P x P rather than critical 

points of dP : TQ — > R relative to P : TQ — > P x P. One can say that we do not work directly 
with TQ but we do it via the map £ 1— > 0^7. 



17 



If T 7 Q and E ro ^ are isomorphic (as we assumed in this paragraph), then the above setting 
coincides with the setting considered in the previous paragraph. However, the results and con- 
structions of this paragraph hold even if the two spaces are not isomorphic. The reason is that 
what one really needs is the image of ker P in TQ under £ i— > ^7 and the image of ker TP in 
TTQ under (A£, 77) H> <5(a£,»j)7- 

Finally let us comment on the isomorphism between E ro ^ and T 7 Q. It was proved in [1] 
that, if E is an AL algebroid, then given an element S^fit) € T^Q and £0 € i? TO7 (f ) such that 
P(£o) = Tr(<57(io)) there exists a unique £(i) E such that ^7 = £7 and £(to) = £0 (to find 
such a £ one needs to solve a linear ODE). If E is not transitive (i.e. ker p is nontrivial) there 
may be several elements £0 and several corresponding curves £(t) for which £7 = ^7. In such a 
case £ 1 — y ^7 is just an epimorphism. 

Symmetry of the second variation. The question if 5 2 Sl('j) is a symmetric bilinear form 
needs some attention. First of all to speak about symmetry we have to restrict our attention to 
pairs of generators vanishing at the end-points, i.e. consider 5 2 S[ J ('y)(-, •) : E® ^ x E® — > R. 

Recall that the second variation 5 2 5*1,(7) (^ £) i s defined by means of a vector #(A£,7y)7 with 
generators A£ and 77 vanishing at the end-points. To define o~ 2 Sl(j)(£,, f]) we need another vector 
^(A?j,£)7 with A77 and £ vanishing at the end-points. For a given (A£, 77) there is a natural choice 
of A77 - we take a unique vector in T V E K-related to A£. The comparison of vectors <5(A£,r/)7 
and 5(Ar;,^)7 can answer the question of symmetry. The conclusion is rather surprising: we can 
guarantee the symmetry of the second variation only for Lie algebroids. Simple examples form 
Section 16.11 show that this condition cannot be omitted. 

Lemma 4.6. Let 7 be the solution of the EL equation (|3.2p on an AL algebroid E. If E is a Lie 
algebroid then the second variation 5 2 Sl{^) is symmetric after restricting to E® Qry x E^. . 

Proof. Take any £, rj € -Ey and A£ € T^E vanishing at the end-points. Denote by A77 € T V E 
a unique vector field along r\ which is K-related to A£. Observe that since k is a linear relation, 
then An also vanishes at the end-points. Now 

5 2 S L ( 7 )( V ,£) = f 1 d T 2 L(S {A ^)dt 

and ^ 

<5 2 Sl(7)(^) = / 1 d T 2 L(5 [Av ^)dt. 

Since d^L is K^-invariant the later equals J^ 1 d"[ 2 L{KE <5(Ar/,§)7)di- Our goal is to compare 
vectors <5(A£,?7)7 and ke <5(A??,£)7 = ^At/OV?)- They both belong to Tg^TE and both project to 
(L7 under Tr^;. We claim that if E is a Lie algebroid then 

<5(A^)7 = <Wo7- (4.4) 

Above equality has a simple coordinate proof. The reader is invited to do the necessary calcu- 
lation using formula (|5. 11)) form the next section. A short calculation shows that the difference 
^(A£,r;)7 — 3(Ari,(,)7 contains the RHSs of equations (|5.ip and (|5.2p which are local versions of 
(|2,2p and (|2.3p . To get the equality these conditions have to be satisfied. In fact this difference 
is a vector in TTE vertical w.r.t. Tte which corresponds to (^7, J{j, Vi £)) i n the canonical 
identification YTE ~ TE x e TE. Here «/(•,-,•) denotes the Jacobiator. 

The geometric interpretation of the above calculations is the following. According to [5] E 
carries the structure of a Lie algebroid if and only if the associated DVB morphism s : T*E — > 
TE* maps the canonical Poisson structure on T*E to the bivector field diTI associated with the 
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algebroid structure d^E, i.e. e is a morphism of algebroid structures. This condition can be also 
expressed in the language of relations k. Roughly speaking Tk should intertwine relation d^K 
and ke- This is the sense of equation (|4.4|) . Some more work is needed to make this intuition 
rigorous. □ 

JVFs and the second variations. Now we are ready to formulate the main result of this 
section. 

Theorem 4.7. Let L : E — > R be a Lagrangian system on an AL algebroid r : E — > M and let 
7 be a solution of the EL equation (|3.2|) . A curve rj 6 E to ^ is a JVF along 7 if and only if 

<5 2 S L ( 7 )(r?,£) =0 for every £ e E^. 

If E is a Lie algebroid then the null space of the symmetric bilinear form 

5 2 S L ( 7 )(;-):E TOy xE TO1 ^R 

consists of all JVFs along 7 vanishing at the end-points to and t\ (i.e. 7(^0) an d l(ti) are 
conjugated alongj if the null space is non-trivial). The multiplicity 0/7^0) and"f(ti) as conjugate 
points equals the dimension of this null space. 

Proof. The proof is straightforward. We know that r\ is a JVF if and only if ^7 satisfies the 
EL-equation for the Lagrangian d^L : TE — > K on the lifted algebroid d^E (cf. Theorem 13 . 6[) . 
From Remark 13.21 the later is equivalent to the vanishing of 

drSd T L(o~Az( S v r Y)) = dTdTS L (5^(S v j)) = d T 2 S L (K E o 6^(6^)) = dT 2 £x((<!>(A£,r;)7)) 

for every d'p-E'-admissible variation 5/±£ (^7) with vanishing end-points. From (|4.2p we know that 
this is equivalent to the vanishing of 5 2 Si J ('y)(r], £) for every £ € E® Q ^. 

The second part of the assertion is trivial. From Lemma 14.61 we know that S 2 Sl^)^, •) is 
symmetric for Lie algebroids and that, by definition, existence and multiplicity of conjugated 
points along 7 is defined by means of JVFs along 7 vanishing at the end-points. □ 

Observe that in general we cannot expect the dimension of the null space of 5 2 Sl('^) to be 
finite (as is in the Riemannian case). Note namely, that Jacobi equation (|3.5p is an implicit linear 
differential equation, so we cannot estimate the number of independent solutions. In particular 
if L = 0, then ever curve £ 6 E rol is a JVF. However under extra assumptions for L it is possible 
to control the number of solutions and consequently, the dimension of the null space. 

5 Coordinate calculations 

So far we worked in purely geometrical terms. In this section we give local description of the 
most important results and formulas. 

Manifolds and vector bundles We will be using the same notation conventions as before; 
i.e. M, Q are manifolds, r : E —> M is a vector bundle ir : E* — > M its dual, etc. 

On Q we introduce local coordinates (q a ) a = 1, . . . , m. We have induced coordinates (q a , qP) 
on TQ and (q a ,P/3) on T*Q. On the iterated tangent bundle TTQ we have induced coordinates 
(q a ,qP,6ql,6q s ); {q a ,pp,q 1 M on TT*Q and (q a , qP, Pj , -k s ) on T*TQ. 

On r : E — >• M we introduce local coordinates (x a ,y 1 ), a = 1, . . . ,n; i = l,...,k, where 
(x a ) are local base coordinates and (y l ) are linear coordinates on fibers. The choice of such 
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coordinates is equivalent to the choice of a basis of local sections (ei, . . . , e k ) of r : E — > M. On 
the dual bundle ir : E* —¥ M we have dual coordinates (x a ,Zi). 

After choosing the basis (ei, . . . , e&) we can introduce natural local coordinates on tangent 
and cotangent bundles of E and E* 

(x a , y\ x b , yi) in TE, (x a , z i: x\ Zj ) in TE*, 

(sV.P&.tt,-) in T*E, (x a ,^,p b ^) in T*E*. 

Algebroid structure. The skew-algebroid structure on E is defined by means of functions 
Pi(x) and c l - k {x) given by 

and 

[ej,ej] = ei ■ Cj k (x). 

The skew-symmetry of [•, •] reads as c l - k {x) = — c^-(x); condition (|2.2p for E being AL takes the 
form 

p a k {x)c%{x) - {^(x)p^) - ^{x)p){x)\ = (5.1) 

ant condition (|2.3[) for E being Lie is (|5.ip together with 

dc s dc s dr s 

c f fe (x) c <.(x)+4(x^^ (5.2) 

The canonical bivector IT £ X 2 (E*) associated with the algebroid structure on r reads as 

1 
2 

and the induced DVB map e : T*E — > TE* as 



U(x, z) = ^4(x)Z k d Zi A d z . + p b (x)d Xb A d Zi (5.3) 



e{x a ,y l , Pb ,Zj) = (x a ,z i ,p b k (x)y k ,c k j (x)y i z k + pj(x)p a ). (5.4) 
Relation k dual to e (see f)2.5|) ) relates vectors 

« : (x a , Y\ p\{x)y\ Y*) - (>,</, o\{x)Y\ V + c^x)/^) . (5.5) 

It follows that if j(t) ~ (x a (t), y l (t)) is an admissible curve and £(t) ~ (x a (t), f l (t)) a generator 
of the infinitesimal variation Sgj (see Definition 12. 3p then 

8a-(x%t),y\t),pt(x(t))f\t)J^t) + 4 l (x(t))y k (t)f l (t)) , 

where j(t) satisfies additionally the admissibility condition 

x a (t)=pf(x(t))y i (t). 

Canonical DVB isomorphisms. For a particular example of the tangent bundle Lie algebroid 
tq : TQ — > Q the coefficients of the anchor and the bracket in canonical coordinates (q a ,q^) 
are simply Pp(q) = 5^ and c^(q) = 0. Hence we have the following coordinate formulas for 
canonical DVB morphism a Q : TT*Q -> T*TQ and kq : TTQ ->• TTQ: 



otQ ■ (q a ,Pi3,f,Ps) ' — > (q a ,q P ,P~,,P6^J , 
k q : (q a ,qP,5q\5q 8 ) .— > (q a , Sq? , q\ Sq° ) . 
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Tangent lifts. For a function F : Q — > R its tangent lift d^F : TQ — >• R reads locally as 

d T f(ff,4) = ^(9)9° 



and d T 2 -F : TTQ ^ R as 



r)F Ffi F 

dT 2 -P(g, q, Sq, Sq) = _^ + ^_(,)^. 



The Kg-invariance of di 2 ^ follows from the symmetry of the second derivative Q~^§^ff(q) = 
S 2 F / \ 

The tangent lift d T n € Af 2 (Ti?*) (see (I2TTT]) ) of the bivector n G ;f 2 (£*) given by (JO 



reads as (see (j2.8j) ) 

1 / 9c fc - \ 

d T n(x, z,x,z) = -[ 4(x)i fe + -^{x)z k x a 8% A <9^ + c^(x)z k d Zi A <9 % + 

/ (5.6) 



EL and Jacobi equations. For a Lagrangian function L : E — >• R the maps and read 



a.s 



Al(x,J/) ~ ^x a ,^(x,y),p b l (x)y k ,c^ j (x)y i ^{x,y) + p^{x)-^{x,y) 

Now the EL equation (|3.2p for 7(2) ~ (x a (t) , y l (t)) takes the following form 

^(t) = pf(x(t))y'(t) (5.7) 
1 (|^(^),y(t))) =^(x(*))|^(x(t),y(*)) + 4(x(<)y(t)||( af (t) > y(*)). (5.8) 



We see that (|5.7p implies admissibility of 7(i). 

The Jacobi equation ([33]) for a JVF 77 (i) ~ (x a (i),f (t)) along 7(4) ~ (x a (t), takes the 

following form 

pt(x)c^(x)y^e + ^Wpfctf? = ^BPfcWe (5-9) 

and 

d ( d 2 L „, d 2 L ( .„ „ , N j , 



sUv' ?( * + wt f+4( * { 



a , , A ( d 4j, x s dL k d 2 L dp) 8L b , , d 2 L \ 

fr _ , 3 2 L _ <9 2 L 



(r + cUx)y i e) {c%{x)^- k + c%{x)y^ k + $( x ) 



QySQ x a 



where j(t) satisfies (|5.7p and (|5.8p . observe that (|5.9p is just (|5.ip multiplied by y 3 and £ l . This 
is in agreement with our considerations form Remark 13.51 
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The second variation. At the end of our considerations we will give the coordinate description 
of vectors <5(A£,r/)7 which generate second variations. We will be using induced coordinates 

(x a ,y\x b ,y^5x c ,Sy k : 5x d ,dy l ) 

on TTE. Our ingredients are: 

• an admissible curve 7(f) ~ (x a (t) , y l (t)) € E i.e. x a = p?(x)y l , 

• curves £(t) ~ (x a (t) , f (t)) and r){t) ~ (x a (t), ti(t)) from E T0J and 

• a curve Af(t) ~ {x a (t), f(t), Ax b (t), Af (t)) 6 T ? £ such that Tr o A( = i.e. Ax a = 
pf(x)h l . 

Now <5(Af,fj)7 equals 

x°, j,*, ^(*)f , /i + ^(x)//, ^(x)^, A fc + c^s)^*, M.p-( x )Vf + pf(x)Af, Ay 1 

(5.11) 
where 

dc l ■ 

Ay 1 := A/' + r| ; i,-),/A/< + + ^(.r^// )/' + ^(x)p%(x)h k y i f . 



6 Examples 

6.1 Riemannian geometry 

E'-connections. Consider an AL algebroid structure on r : E — > M with the anchor p and the 
bracket [•,•]. An E -connection is a bilinear map 

V : Sec(£) x Sec(£) — ► Sec(£) 

satisfying the following properties: 

V fX Y = fV x Y 

V x fY = fV x Y + p(x)(f)Y, 

for every Sec(-E) and / £ C°°(M). Observe that even though V is defined on sections of 

E, to compute V x Y&t p € M it is enough to know X(p) and the derivative of Y in the direction 
P(X)(p). 

An ^-connection is called torsion-free (compatible with the algebroid structure) if 

V X Y- V Y X = [X,Y]. 

If E is equipped with a Riemannian metric, i.e. a fiber-wise linear symmetric and non-degenerate 
map g : E Xm E — > M., we can speak about metric ^-connections which satisfy 

p(X)g(Y, Z) = g(V x Y, Z) + g(Y, V X Z). 

The following well-known result generalizes the fundamental fact form Riemannian geometry. 
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Lemma 6.1. Let (E,g) be as above. Then the formula: 

2g(V x Y, Z) :=p(X)g(Y, Z) + p(Y)g(X, Z) - p(Z)g(X, Y)+ 
g([X, Y],Z)- g([Y, Z],X)+ g([Z, X],Y) 

defines a unique metric and torsion-free E -connection. 
Given an S-connection V we can define its curvature: 

R V (X,Y)Z := V x VyZ- V Y V X Z - V [xy] Z. 

A straightforward computation shows that R\?(X, Y)Z is tensorial in the first two entries and 
tensorial w.r.t. Z whenever E is an AL algebroid. 

A geodesic problem. Given E and g as above we can consider a natural Lagrangian system 
L : E R on E given by the energy function L(a) := \g(a, a). Our goal is to describe the EL 
and the Jacobi equation for this system. We will show that the EL equation for an admissible 
curve 7 G E reads as 

V 7 7 = 

and the Jacobi equation for rj £ E TO ~ f as 

V 7 V 7 7 ? + J R(r ? , 7 ) 7 = 0. 

Of course in the above equations V is the canonical connection described in Lemma 16.11 and R 
its curvature, obviously these equations are straightforward generalizations of classical formulas 
from Riemannian geometry. 

Calculations. In our considerations we will be working with geometric objects defined along 
an admissible curve "f(t) G E. However sometimes we will be using objects defined locally around 
7(t), such as the algebroid bracket [•, •] and the covariant derivative V. Therefore we will need 
the following construction: if £(t) is a curve in E over x(t) G M and A£ £ T^E a vector field 
along £, then by a ^-extension of £ we will understand any local section £ : M — > E around £ 
such that l(x{t)) = {{t) and T£ (Tr(A£(i))) = Af (t). 

In particular let "f(t) G E be an admissible curve over x(t) G M and r\ G E to1 a generator of 
an admissible variation ^7. Then for 7 - a (^7-extension of 7 and r] - any extension of 77 we got 

D?>7] = along x(t). (6-1) 

This fact is just a matter of a simple coordinate calculation. 

Now we are ready to derive the EL equation. We claim that for an admissible curve -y(t) G E 
we have 

XMt)) =g( 7 (t)) (6.2) 

and 

Aj(7(*)) = T?(<y(t) - F( 7 (t), V 77 (<))) , (6.3) 

where g : E 3 a4 g(a, ■) G E* is a natural isomorphism and V(a, b) G X a E C T a E is a vertical 
vector corresponding to (a, b) in the canonical identification YE ~ E x^/ £\ Note that (|6.3p is 
tensorial w.r.t. 7 although the RHS contains 7. 
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Formula (|6.2p can be checked directly in coordinates. To prove (|6.3p we will do the following 
reasoning. For any admissible j(t) € E over x(t) G M and any generator G E TQ ^^ t y. 



1 



(dS L (j),6tf)= (dL( 7 (t)),5 a (t))dt = -Tg(6 a (t),5 a (t))dt = -p(d)g(j,j)(x(t))dt, 
where 7 is an ^7-extension of 7. The later equals 

<7(Vf7,7)(*(*))d< = 7] + V 7 e,7)(x(t))d< ^ T <7(V 7 e,7)(*(*))d* = 

to J to ■/ to 

1 p(7)5(e,7)(^(*))di - / 1 g(£, V 77 )(x(i))di. 



'to J to 

Since 7 is admissible p(j)g(£, 7 )(x(i)) = g7<7(£(i),7(i)); and hence 



<d5 (7 ),^7> =<?(£(*), 7(*)) 



ti /-ti 



5 (e,v 77 )(x(t))dt 



<eW,A L (7(y))) 



to J ta 



s(£,V 77 )(s(t))dt. 



to ■ / *o 

Comparing the above with (|3.3p we get 

(e(*).V AL(7(t)) _ d AL(7(t)) > = -gift), V 7 7(*))- 

Obviously 

<^)<Hv L ( 7 (t))-A Ai(7{T)) } = -</(£(*), V 77 (t)) = (e(*),?(-V 7 7(t))) = 

(£(<), V Tff( 7 (t)-F( 7 (t),V 77 (t))- 7 (t))) = (^(*)' V Tg( 7 (t)-V( 7 (t),V 77 (t)))-^A £ ( 7 (t)))- 
It follows that (|6.3p holds. In particular we know that the EL equation reads as 

Tflf(^( 7 (t),V 77 (t)))=0, 

which is equivalent to V 77 (i) = 0. 

Now we are ready to derive the Jacobi equation. According to (|3.5p we need only to compute 
the tangent maps of Ax and Al at <L 7 . Obviously 

TA L (A 7 ) =Tgf(A 7 ). 

Now 

k e * o TA L (A 7 ) =k e * o TTg (A( 7 ) - AV(-y, A t7 )) = 

TTg o ke (A( 7 ) - AF (7, A 77 )) = TTg ((A 7 )' - V(A 7 , A(Vyy))) • 

We deduce that the Jacobi equation for 77 E £V 7 along a geodesies 7 takes the form TTg (V(A7, A(V 7 7))) : 
at A7 = 5 v j, which is equivalent to 

A(V 77 ) = 0, 

where A7 = ^7 and V 7 7 = 0. Now if A7 = 5 v j then 

(A£, V TW(A7iA( v 77 )))) = Tg{A£, A(V 77 )) = p(rj)g(£, Vyy), 
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where 7 is a ^7-extension of 7 and £ is a A£-extension of £. The above equals 

g(V v t, V 7 7) + g(r], V^V^t") = g{i], V v \7^y), 
since V 7 7 = 0. Using the definition of the curvature we get 

g (£, V^V T 7) = 5 (£, 7)7 + V 7 V^ + V^]7) = 9 (£, R(v, lh + V 7 V^) = 
g (£, flfa, 7)7 + V-V-r/ + V-[r/, 7]) ^ g (£, i?(r?, 7)7 + V 7 V 7 t?) 
In this way we have proved that 

#(V,7)7 + V 7 V 7 r? = 

is equivalent to the Jacobi equation and, simultaneously, we have proved the generalization of 
6.2 Euler-Poincare equations 

Tangent lift dxfl. Let (g, [•, •]) be a skew-symmetric algebra. The vector space g treated as 
a vector bundle over a point, together with the bracket [■,•] and the trivial anchor is an AL 
algebroid (condition (|2.2p is trivially satisfied). It is a Lie algebroid if (g, [•, •]) is a Lie algebra. 
We shall now describe the lifted algebroid structure dxg. It is also given by a skew-symmetric 
algebra structure on Tg, as the tangent lift of the vector bundle g — > {*} is again a VB over 
a point. Using the canonical identification Tg ~ g x g we obtain a skew-symmetric algebra 
structure (g x g, [-, -]d T )- 

The canonical isomorphism e : T*g ~ g x g* — > g* x g* k, Tg* for the algebroid (g, [•, •]) 
reads as 

e:(a,z)^(z,(z,[a,-})); (6.4) 

where {■,•) '■ Q* x g — > IR is the canonical parring. It follows that for an admissible curve 
7 = a(t) € g (since the anchor is trivial every path is g-admissible) and a generator £ = f(t) £ g 
the associated admissible variation reads as 

5a = (a(t)J(t) + [a(t),f(t)]) € g x g « Tg. (6.5) 

Proceeding analogously as in Section 2] we can differentiate a 1-parameter family of such variations 
to get 

A(S a ) = (aj+ [a, /], Aa, A/ + [Aa, /] + [a, A/]) (t) € g x g x g x g « TTg. 

On the other hand we know form Section 2] that A(5^j) is just K g o 5^A^y, where 5a^Aj is a 
dxg-admissible variation along A7 = (a(t), Aa(t)) generated by A£ = (f(t), Af(t)). It turns out 
that 

o A5 A 7 = (a, Aa, / + [a, /], A/ + [Aa, /] + [a, A/]) (t)Ggxgxgxg« TTg. (6.6) 

Formula (|6.5p for (g x g,[-,-]d T ) describes the relation between the skew-algebra bracket and 
admissible variations. Hence from (|6.6p we deduce that 

[(a, Aa), (/, A/)] dT = ([a, /], [Aa, /] + [a, A/]) . (6.7) 
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The EL and Jacobi equations. Let us now describe the EL and Jacobi equations for a 
Lagrangian system L : Q — ¥ R on 0. Observe that dL : g — > T*g ~ g x g* defines a map 
5L : g — > g* given by dL(a) = (a, 5L(a)) (in fact 5L = \l in the terminology of Section[3]). From 
(|6.4p we deduce that the EL equation reads as 

±(5L(a(t))) = (5L(a(t))Mt)r])- (6.8) 

After introducing linear coordinates i = 1, . . . , k on Q (denote by c l - k the coefficients of the 
bracket in these coordinates) this equation takes the form 

We recognize the Euler-Poincare equation. 
Since d^L(a,b) = (<5L(a),&) then 

<5d T L(a, b) = (SL(a), 5 2 L(a, b)) € g* x g* « Tg*, 

where 5 2 L : g x g — > g* is a map linear in the second entry. Formula (|6.8[) together with expression 
(|6.7p will give us the EL equation for the lifted system d^L : Tg ->Ion dTg- Then form Theorem 
13.61 we can deduce the Jacobi equation for a Jacobi vector field (a(t),h(t)) £gxg» Tg along 
a(t) € g: 

±{6L{a)) = (6L(a),[a,-]) 

! (6.9) 
— [d 2 L(a, h + [a, /»])) = (SL(a), [h + [a, h], •]> + (5 2 L(a, h + [a, h}), [a, •]>. 

The first equation is just the EL equation for the Lagrangian system L : g — > R. The second 
equation in coordinates a(i) ~ (y l (t)) and /i(t) ~ (h l (t)) reads as 

d / 8 2 L \ ( k dL . , fc 4 d 2 L \ TS 

where J s = h s + c s kl a k h l . In a particular case = |<?(a, o), where g:gxg— >• R is a symmetric 
2-tensor, we get 5L(a) = g(a, •) and 5 2 L(a, b) = g(b, •). In this situation equation (|6.9p takes the 
form 

g(a, •) = g(a, [a, •]) 

g(h + [a, h) + [d, h], •) = g{a, [h + [a, h], •]) + g(h + [a, h], [a, •]). 

Symmetry of the second variation. Finally let us study the symmetry of the second vari- 
ation of the action ^ 

997(^^(7)= / 1 £(7(*))dt- 
J t 

We know that for i] and £ vanishing at the end-points 

<5 2 L( 7 ) (??,£) = f 1 d T 2 L(5 (A ^ v) )dt, 
where A£ vanishes at the end-points. Since S^a^,t])1 = k q° ^A^(^tjT) from fj6.5[) and ()6.6|) we get 
<5(A 5 ,r,)7 =(«,/+ [a, /U + [a, h],Af + [h+ [a, /»],/] + [a, A/]) (t), 
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where 7 = a(t), 77 = h(t), £ = f(t) and A£ = (/(t), A/(t)). The curve Arj = (h(t),Af(t) + 
[h(t),f(t)]) is K-related to A£ and also vanishes at the end-points, since A£, £ and r] do. We 
have 

W)T = + [°. ^ / + [°. /], A/ + [A, /] + [Zi, /] + [/ + [a, /U] + [a, A/ + [/i, /]]) (t). 
A short calculation shows that 

^(Ae,T?)7 - K s <W?)7 = [a,f],0,J(a,h,f)j (t), 

where J(a, Zi, /) := [[a, /i], /] — [a, [/i, /]] + [h, [a, /]] is the Jacobiator. Now 

S 2 S L (lM®-S 2 SL(m,v) = d T 2 L(5 (A ^ )7 )dt - [ tl d T 2 L(5 (Av ^)dt = 

Jt Jto 



/•tl />tl Ptl 

/ d T 2 L(5 (A ^)7)dt - / d T 2 L(K o 5 (Ar?iC) 7)dt = / (dd T L(5^), <5(AC,r,)7 - « 9 <J(A^)7>di 
r<dL( 7 ),J(7,i7,0>d<- 

Jin 



We see that if (g, [•, ■]) is not a Lie algebra (i.e. the bracket [•, •] does not satisfy the Jacobi 
identity) then for a generic Lagrangian 8 2 Sl{^) will not be symmetric. 

Final remarks 

Further study, Hamiltonian dynamics. Let us note that the classical geodesic equation 
can be regarded not only as a EL equation on TM, but also as a Hamilton equation on T*M. In 
fact most of the results of this paper can be repeated for Hamiltonian systems on algebroids. The 
results of Theorems 11.21 and 11.31 hold in the Hamiltonian case after substituting "EL equation" 
by "Hamilton equation" and changing the standard variational principle to a more suitable one. 
A detailed study of this case will be given in a forthcoming publication [9]. 



Curvature? It would be interesting to extend other aspects of the classical theory of Jacobi 
fields and conjugate points to the geometric setting presented in this paper. In particular it 
is tempting to find an analog of curvature for arbitrary Lagrangian (or Hamiltonian systems) 
and use it to prove results about existence or non-existence of conjugate points (see e.g. [7] 
for such results with applications to general relativity). The results of this kind (especially in 
the Hamiltonian context) could be useful in the study of Hamilton- Jacobi-Bellman equation in 
optimal control theory (see e.g. [1]). It seems to us that similar ideas could be used to study 
phase transitions in thermodynamics. 

So far we were unable to find any satisfactory definition of curvature. It is highly probable 
that a general definition of this kind does not exists and some additional assumptions about 
the nature of the systems in consideration should be made. For example the yet unpublished 
paper [8j gives a definition of the curvature for Hamiltonian systems if the system satisfies certain 
geometric conditions. In the simplest case of a Hamiltonian system on a cotangent bundle T*M 
this condition is the hyper-regularity of the Hamiltonian. 
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